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ABSTRACT 



Aims. We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation 
(Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to Eiso ^ io^^~^^ ergs, 
providing a unique probe of the high-redshift Universe. 

Methods. We consider both the so-called Pop III.l stars (primordial) and Pop III. 2 stars (primordial but affected by radiation 
from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical 
evolution of the intergalactic medium. 

Results. We show that Pop III. 2 GRBs occur more than 100 times more frequently than Pop III.l GRBs, and thus should be 
suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained 
by the current radio transient searches. We expect ~ 10 — 10* radio afterglows above ~ 0.3 mjy on the sky with ~ 1 
year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while 
we expect to observe maximum of < 20 GRBs per year integrated over at 2 > 6 for Pop III. 2 and N < 0.08 per year 
integrated over at 2 > 10 for Pop III.l with EXIST, and TV < 0.2 for Pop III. 2 GRBs per year integrated over at z > 6 with Swift. 
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1. Introduction 

The first stars in the Universe are thought to have played 
a crucial role in the early cosmic evolution by emitting 
the first light and producing the first heavy elements 
(iBromm et al. I2OO9I) . Theoretical studies based on the 
standard cosmological model predict that the first stars 
were formed when the age of the Universe was less than 
a few hundred million years and that they were predom- 
inantly massive (lAljel et al.ll2002l: lOmukai fc Pallall2003t 



Yoshida et alll2006[) . Although there are a few observa- 



tional ways to probe the early stellar populations, there 
has not yet been any direct observations of the so-called 
Population III (Pop III) stars. 

Recently, it has been proposed that massive Pop III 
stars may produce coUapsar gamma-ray bursts (GRBs), 
whose total isotropic energy could be i^iso ^ 
ergs, i .e., > 2 orders o f ma g nitude above av 

I2OIOI: 



Meszaros fc Reedl20ld iBarkovllioiol: ISuwa fc Iokall201ll ). 

Even if the Pop III star has a supergiant hydrogen en- 
velope, the GRB jet can break out the first star be- 
cause of the l onglasting accretion of the envelope itself 
(jSuwa fc Iokall201ll : iNaeakura et"alll201l[ ). The observa- 
tions of the spectral peak energy-peak luminosity rela- 
tion also imply that the GRB for mation rate does not 
significantly dec rease to 2; ~ 12 (jYonetoku et al ] I2OO4I: 
Nava et al.i[2oTlh . 

Observations of such energetic GRBs at very high red- 
shifts will provide a unique prob e of the high-redshift 
Universe (iLamb fc ReichartI l200d iGiardi fc Loebl I2OOOI: 



Gou et al.lT2004t lYonetoku et al.ll2004l ). Iloka fc Meszaros 
(|2005h study the radio afterglows of high-redshift GRBs 
and show that they can be observed out to z ~30 
by VLAB, LOFA RgI and SKAI (sec also |lo3 l2003l) . 
Inoue et al.l (|2007l ) calculate time-dependent, broad-band 



erage (jKomissarov fc Barkov 



Toma et al. 2011 



http://www.vla.nrao.edu/ 
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afterglow spectra of high-redshift GRBs. They suggest 
that spectroscopic measurements of molecular and atomic 
absorption lines due to ambient protostellar gas may be 
possible to z ~ 30 and beyond with ALMAQ, EVLAH, and 
SKA. In the future, it will be also promising to observe the 
GRB afterglows located by gamma-ray satellites such as 
SwifB, SVOMH, JANUS Hand EXISt|1. Clearly, it is im- 
portant to study the rate and the detectability of Pop III 
GRBs at very high redshifts. 

There have been already a few observa tions of GRBs at 
high redshifts. GRB 090429B at z = 9.4 (jCucchiara et al 



20111 ) is the curre nt record-holding ob ject, followed by 
a 2 = 8. 6 galaxy dLehnert et alll2010l). GRB 090423 at 



26 ( Salvaterra et al 



2009; 



Tanvir et alll2009[) . GRB 
200^7gRB 050904 



080913 at z = 6.7 (iGreiner et al _ 

at z = 6.3 (|Kawai et all 120061: iTotani et al.l l2006l). and 
the h ighest redshift quasars at z = 7.085 (IMortlock et al 
1201 ll) and z = 6.41 (jWillott et alll2003l ). 1c'handra et al' 
(I2OIOI ) report the dis covery of radio a fterglow emission 



from GRB 090423 and lFrail et all (|2006l) for GRB 050904. 
Observations of afterglows make it possible to derive the 
physical properties of the explosion and the circumburst 
medium. It is intriguing to search for these different sig- 
natures in the GRB afterglows at low and high redshifts. 

The purpose of the present paper is to calculate the 
Pop III GRB rate detectabl e by the current an d fu- 
ture GRB missions (see also ICampisi et al.l I2OIII) . We 
con sider high-redshift G RBs of two populations follow- 
ing 



Bromm et al 



()2009l) . Pop III.l stars are the first- 
generation stars that form from initial conditions deter- 
mined cosmologically. Pop III. 2 stars are zero-metallicity 
stars but formed from a primordial gas that was influenced 
by earlier generation of stars. Typic ally, Pop III. 2 stars 
are formed in an initially ionized gas (jjohnson fc Bromm 
l2006l:lYoshida et al.ll2007h . The Pop III.2 stars are thought 
to be less massive (~ 4O-6OM0) than Pop III.l stars (~ 
IOOOMq) but still massive enough for producing GRBs. 

We have calculated the GRB rate for these two popu- 
lations separately for the first time. The rest of the paper 
is organized as follows. In Sect. 2, we describe a semi- 
analytical model to calculate the formation rate of pri- 
mordial GRBs. In Sect. 3, we show our model predic- 
tions and calculate the detectability of Pop III GRBs 
by future satellite missions and by radio observations. In 
Sect. 4, we discuss the results and give our concluding 
remarks. Throughout the paper wc adopt the standard 
A cold dark matte r model with the b est fit cosmologi- 
cal parameters from lJarosik et al. I (l201lh (WMAP-YrTEB), 
fim = 0.267, nA = 0.734, and Hq = 71km s'^Mpc^^ 



* www.alma.nrao.edu/ 

^ http://www.aoc.nrao.edu/evla/ 



2. Gamma-ray burst rate 

We assume that the formation rate of GRBs is propor- 
tional to the star formation rate ( Totanjri997 : Ishida et al 



20111 ) . The number of observable GRBs per comoving vol- 
ume per time is expressed as 



*grb(^) 



4tt 



-?7GRB Vhcam^*{z) 



logilim(2) 



p{L)d\ogL, 

a) 

where t?grb is the GRB formation efficiency (see section 
2.6), T^beam the beaming factor of the burst, flots the field 
of view of the experiment, ^» the cosmic star formation 
rate (SFR) density, and p{L) the GRB luminosity function 
in X-rays to gamma rays. The intrinsic GRB rate is given 

by 



GRB 



(z) = ??GRB**(^)- 



(2) 



The quantity iiim(z) is the minimum luminosity threshold 
to be detected, which is specified for a given experiment. 
The nonisotropic nature of GRBs gives Tyboam ~ 0.01 — 
0.02 (iGuetta et aD 2005 ). Using a radio transient survey. 
Gal- Yam et al.l ( 20061 ) place an upper limit of ?7bcam ^ 
0.016. Given the averag e value of jet opening angle 9^6° 



(jGhirlanda et al.l |2007^ and r^b. 



1 5.5 X 10 we set 
Vheam = 0.006 as a fiducial value. The adopted values of 
17obs are 1.4, 2, 4, a nd 5 for Swift, SVOM JANUS, and 
EXIST, respectively (ISalvaterra et al.ll2008l ). 



2.1. The number of collapsed objects 

We first calculate the star formation rate (SFR) at 
early epochs. Assuming that stars are formed in col- 
lapsed dark matter halos, we follow a popular prescrip- 
tion in which the number of c ollapsed objects is calculated 
by the halo mass function dHernquist &: Soringell 12003 : 
Greif fc Bro"^l2006l: iTrenti fc Stiavellill2009l). We adopt 
the S heth-Tormen mass function, /st, (jSheth fc Tormen 
19991) to estimate the number of dark matter halos, 
nsT{M, z), with mass less than M per comoving volume 
at a given redshift: 



fsT = A 



■ exp 



' 2a2 



(3) 



where A = 0.3222, ai = 0.707,p = 0.3 and 6c = 1.686. 
The mass function /st can be related to the nsT{M, z) as 



/sT = 



M dnsT{M,z) 



Pn 



dlnc7-i 



(4) 



where pm is the total mass density of the background 
Universe. The variance of the linear density field (j{M, z) 
is given by 



http;/ /swift. gsfc.nasa.gov/docs/swift/swiftsc. html 
^ http://www.svom.fr/svom.html 

* http:/ /sms. msfc.nasa.gov/xenia/pdf/CCE2010/Burrows.pdf 
^ http:/ /exist. gsfc.nasa.gov/design/ 

http:/ /lambda.gsfc. nasa.gov/product/map/current/ 



cr2(M,z) 



27r2 



k^P{k)W^{k,M)dk, (5) 



where b{z) is the growth factor of linear perturbations 
normalized to 6 = 1 at the present epoch, and W{k, M) 
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is the Fourier-space top hat fiher. To calculate the power 
spectrum P{k), we use the CAMB codcEl 

for our assumed 

ACDM cosmology. 

2.2. H2 Photodissociation 

The star formation efficiency in the early Universe largely 
depends on the ability of a primordial gas to cool and con- 
dense. Hydrogen molecules (H2) are the primary coolant 
in a gas in small mass "minihalos" , and are also fragile 
to soft ultraviolet radiation, and thus a ultraviolet back- 
ground in the Lyman- Werner (LW) bands can easily sup- 
press the star formation inside minihalos. We model the 
dissociation effect by setting the mi nimum mass for hal os 
that arc able to host Pop III stars (lYoshida et al.ll2003[ ). 

For the minimum halo mass capable of cooling by 
molecular hydrogen in the presence of a Lyman- Werner 
(LW) background, we adop t a fitting for r nula g iven by 



Machacek et al.l ( 2001 ) and 



Wise fc Abej (120051). wh ich 



also agrees with results from 



O'Shea fc NormanI (|2008h : 



exp 



fed 

0.06 



(1.25 X 10^ 



.7xlO'F^4%,), (6) 



where i^LW,-2i = 47rJLw is the flux in the LW band in 
units of 10^^^erg^^s~-'^cm~^Hz~^, and fed the fraction of 
gas that is cold and dense. We set fed = 0.02 as a con- 
servative estimate. We compute the LW flux consistently 
with the como ying density in stars Pi .(z) via a conversion 
efficiency ryLw (| Greif fc Bromm l2006h : 



he 



(7) 



Here, ?7lw is the number of photons emitted in the LW 
bands per stellar baryon and toh is the mass of hydrogen. 
The value of t^lw depends on the characteristic mass of 
the formed primordial stars, but the variation is not very 
large for stars with masses greater than ten solar masses 
(|Schaerej|2(303) . We set 77LW = 10'' for both Pop III.l and 
Pop III. 2 for simplicity. 

Next we calculate the stellar mass density as 



dt 



dz' 



dz'. 



(8) 



For a given z, the integral is performed over the maxi- 
mum distance that an LW photon can travel before it is 
redshifted out of the LW bands. The mean free path of 
LW photons at z = 30 is ^ 10 Mpc (phys ical). Photons 
trave l over the mean free path in ^ 10^ yr (jMackev et al 



2003h . Halos with virial temperature less than 10* Kelvin 
cool almost exclusively by H2 line cooling, and produce 
mostly massive stars. We adopt the mass of such halos, 
M(T'vir = lO^K, z), as an upper limit of halos that pro- 
duce Pop III.l stars. In larger halos, the gas is ionized at 
virialization, and thus the formed stars have, according 
to our definition, similar properties to Pop HI. 2 stars. We 



http:/ /camb.info/ 



assume that A/(Tvir ~ lO^K, z) is the minimum halo mass 
for Pop HI. 2 star formation. 

The collapsed fraction of mass, Fco\{z), available for 
Pop HI star formation is given by 



col W = — / dMMnsT{M,z) (9) 

Pm J Ma, 



piu.i 



for Pop III.l stars, and 



1 



coi(^) = — / dMMnsT{M,z) (10) 

for Pop HI. 2. Using the above criteria, the SFR of Pop HI 
stars can be written as 



III.l 



(z) = (1 - Qmi{z))il - C{z, V^ind))Pmfbf>. 



dF 



III.l 

col 



dt 



for Pop III.l stars, and 



III. 2 



(z) = Qmi{z){l - Ciz,Vy,ind))Pmfbf, 



jpIII.2 

col 

dt 



(11) 



(12) 



for Pop HI. 2. Here, Ci^^'^wind) represents the global fill- 
ing fraction of metals via galactic winds (see section 2.4), 
Qhii(-z) the volume filling fraction of ionized regions (see 
section 2.3), and ft is the baryonic mass fraction. For the 
star formation efficie ncy, we use the value = 0.001 as a 
conservative choice ([Greif fc Broimn 2006) and /* =0.1 
(jBromm fc Loebl200a as an upper limit. The latter choice 
is not strictly consistent with the assumption made in Eq. 

We explore a model with = 0.1 simply to show a 
very optimistic case. 

2.3. Reionization 

Inside growing Hii regions, the gas is highly ionized, and 
the temperature is ~ 10"* K, so the formation of Pop III.l 
stars is terminated according to our definition. The for- 
mation rate of Pop III.l is reduced by a factor given by 
the vo lume filling fraction of ionized regions, Qmiiz). We 
follow IWvithe fc Loeb ( 20031 ) to calculate the evolution of 
Qhii(^:) as 



dQ 



HII 



A^ion dR 



col 



dz 

whose solution is 



0.76 dt 



^ ,n 

a-' 



where 



and 



F{z',z) = - 



2 obu^ 



■C[f{z')-fiz)l 



/(z) = W(l + z)3 + 



1-0. 



(13) 



(14) 



(15) 



(16) 



Here we have assumed the primordial fraction of hy- 
drogen of 0.76. In the above equations, = N^f*fosc is 
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Then we can express fchem, the ratio of gas mass enriched 
by the wind to the total gas mass in each halo, as 



/chcm(A'/, 2, Wwind) = -TT ■ 
6 Vu 



(18) 



where Vh oc is the volume of each halo. The halo radius 
R}i can be approximated by 



UI 



Att X 180p„ 



1/3 



(19) 



Fig. 1. Reionization history calculated using our model. 
The blue line is our model prediction, and the dotted black 
line is the best fit of CAMB code. 



an efficiency parameter that gives the number of ionizing 
photons per baryon, where f^sc is the fraction of ioniz- 
ing photons able to escape the host galaxy, and is the 
time averaged number of ionizing photons emitted per unit 
stellar mass formed. The quantity = 1.95 x 10~^cm~^ 
is the present-day comoving number density of hydrogen, 
ae = 2.6 x 10~^^cm^s~^ is the hydrogen recombination 
rate, and C = (^h)/"-h the clumping factor . We use the 
average value C = 4 (see Pawlik et al.l (|2009f ) for detailed 
discussion about redshift dependence of C). We set the 



Equation ([T5)) takes the self-enrichment of each halo into 
account. The next step is to evaluate the average metal- 
licity over cosmic scales. The fraction of cosmic volume 
enriched by the winds can then be written as 



Ciz,v- 



wind / — 

Pin 



dMf],f^fchcm{M, Z, 



values Z,.,^ = 0.7, /h. = 0.01, a nd iV^ = 9 x 10"^ as fiducial cal metalhcity (see e.g. [Schneider et al 



ad)MnsT{M,z). 

(20) 

Although this may appear a significant oversimplification, 
the model with Wwind as a single parameter indeed provides 
good insight into the impact of the metal enrichment. 

We adopt three different values of Vvrind and exam- 
ine the effect of metal enrichment quantitatively. We as- 
sume that Pop III stars are not formed in a metal-enriched 
region, regardless of the actual metalhcity. Even a sin- 
gle pair instability supernova can enrich the gas within 
a small halo to a met alhcity level we l l abov e the criti- 

2QQ(i)). We effec- 



values (|Greif fc Brommll2006l) . 

In Fig. [T] we show the reionization history calculated 
using our model in comparison with a fitting function that 
is the default par ametrization of reionization in CAMB 
(|Lewis et al.ll2000t) . 



tivel y assume that the so-called cr i tical metalhcity is very 



low (Schneider et al. 


2002, 2003 


Omukai et al. 2005: 


Frebel et al. 


2010). 



20071: iBelczvnski et al. 



2.4. Metal enrichment 

We need to consider metal enrichment in the intcrgalactic 
medium (IGM) in order to determine when the formation 
of primordial stars is terminated (locally) and when the 
star formation switches from the Pop III mode to a more 
conventional one. 

It is thought that Pop III stars do not generate strong 
stellar winds, and thus the main contribution to the 
metal pollution com es from their supernova explosions. 
Madau et al.l (j2001l ) argue that pregalactic outflows from 
the same primordial halos that reionize the IGM could 
also pollute it with a substantial amount of heavy ele- 
ments. To incorporate the effect of metal enric hment by 
galactic wi nds, we adopt a simfiar pre scription to I Johnson 
(|2010[) and lFurlanetto fc Loebl (|2005[) . 

We assume that star-forming halos ( "galaxies" ) launch 
a wind of metal-enriched gas at ^ 20. The metal- 
enriched wind propagates outward from a central galaxy 
with a velocity Wwind, traveling over a comoving distance 
-Rwind given by 



Figures and [3] show the star formation rate (SFR) 
history for both Pop III.l and Pop III. 2 considering three 
different values of the galactic wind, v^md = 50, 75, 100 
km/s. Figure [5] shows that the metal enrichment has lit- 
tle influence on Pop III.l. This is because Pop III.l for- 
mation is terminated early due reionization. In Fig. [3] 
we compare the Pop III. 2 SFR hist ory with a compila- 



tion o f independent measures from iHopkins fc Beacom 
(|2006l ) up to z « 6 and from observations of color- 



selected Lyman bre a k gal axies ( Mannucci et al.l 2007 



Bouwens et al.l l2008l 120111) . L yg Emitters (lOta et al 



20081 ). UV-+TR me asurements ( Reddy et al.l l2008l) and 
GRB observations (jCharv et al.ll2007l: lYfiksel et al.ll2008 : 
Wang fc Dai 2009[) at higher z (hereafter, these will be re- 
ferred to as H2006, M2007, B2008, B2011, O2008, R2008, 
C2007, Y2008, and W2009, respectively). The optimistic 
case for Pop III. 2 is chosen to keep the SFR always below 
the observationally determined SFR at z < 8. 

We compared our model results wit h the SFRs esti- 
mated by other authors in the hterature Bromm fc Loebl 
(|2006l) . lTornatore et al.l (l2007l). iTrenti fc Stiavellil (|2009[ ). 
and also see Naoz fc Bromberg ( 2007t) . It is important 
to note that Pop III formation can continue to low red- 
shifts (z < 10) depend i ng on the level of metal enrich- 



R 



dt 



■wind — / '^wim 



.{l + z')-dz' 



(17) 



ment. iTornatore et al.l ([20071) use cosmological simula- 
tions to show that, because of limited efficiency of heavy 
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Fig. 2. Pop III.l star formation rate. Calculated for for 
weak and strong chemical feedback models and a moder- 
ate star formation efficiency with /* = 0.05. The results 
are shown for Wwind = 50 km/s , red line; and 100 km/s, 
blue line. We also show the theoretical SFRs in the liter- 
ature, from Bromm & Loeb 2006 (Pop Ill.l-f III.2), dot- 
ted black line; Trenti & Stiavelli 2009 (Pop III.2), dashed 
orange line; and Tornatore et al. 2007 (Pop III.l-hIII.2), 
dot-dashed brown line. The purple line is our optimistic 
model where we assume a very high star formation effi- 
ciency, ^ 0.1, and low chemical enrichment, fwind — 50 
km/s. 




5 10 15 20 



z 

Fig. 3. Pop III. 2 star formation rate. Calculated for three 
different chemical feedback models; i^wind = 50 km/s, red 
line; Wwind = 75 km/s, blue line; and tiwind = 100 km/s, 
green line. We also show the theoretical SFRs in the liter- 
ature, from Bromm & Loeb 2006 (Pop III.l-hIII.2), dot- 
ted black line; Trenti & Stiavelh 2009 (Pop III.2), dashed 
orange hne; and Tornatore et al. 2007 (Pop III.1-I-III.2), 
dot-dashed brown line. The purple line is our optimistic 
model where we assume a very high star formation effi- 
ciency, ^ 0.01, and low chemical enrichment, fwind — 
50 km/s. The light points are independent SFR determi- 
nations compiled from the literature. 



element transport by outflows. Pop III star formation 
continues to form down to z = 2.5 (which intriguingly 
matches our model with z^yyi nd = 50 km/s in Fig. [3]). 



The SFR of iTornatore et all (|2007t) has a peak value of 
lO-^Mgyr-^Mpc"^ at z « 6. 

In Fig. U wc also show the result of our model with 
= 0.1 - 0.01 and fwind = 50 km/s for both Pop III.l 
and Pop III. 2, respectively. This model provides an "opti- 
mistic" estimate for the detectable GRB rate for the future 
missions (see Sect. 3). 




2.5. Luminosity function 

The number of GRBs detectable by any given instrument 
depends on the instrument-specific fiux sensitivity thresh- 
old and also on the intrinsic isotropic luminosity function 
of GRBs. For t he latter, we adopt the power -law distribu- 
tion function of Wanderman fc Piran (l2010h 



p{L) = 



-0 2+°-^ 



L < 



,+0.3 



(21) 



L> 



where is the characteristic isotropic luminosity. We set 
L* ~ 10^^ ergs/s for Pop III.l, whereas L* ^ IQ^^ erg sis 
for P op III. 2 stars are similar to ordinary GRBs ([Lil 



2008t I Wanderman fc PiranI [2OI0I) . The Pop III.l GRBs 



are assumed to be energetic with isotropic kinetic energy 
E,so 1056-57g^^ i^yi- long-lived T90 - 1000 s, so that the 



luminosity would be moderate i, ^ e-y x 10^^~^^/1000 ^ 
2^052-53 cj-gg^g if £^ ^ 0.1 is the conve rsion efficiency from 
the jet kinetic energy to gamma rays (jSuwa fc Iokall201l[ ). 



Fig. 4. Comparison of the star formation rates for Pop 

111.1 (blue dotted line) and for PopIII.2 (dashed black 
line), for our optimistic model with a high star forma- 
tion efficiency /* = 0.1 for Pop III.l, = 0.01 for Pop 

111. 2 and slow chemical enrichment Wwind = 50km/s. 

Using the above relation we can predict the observ- 
able GRB rate for the Swift, SVOM, JANUS, and EXIST 
missions. For Swift, we set a bolomet ric en ergy flux limit 



-8 



erg cm 



i^iim = 1.2 X 10 
similar limit for SVOM 
Fiim ~ lO^^erg cm~^ s~ 
nosity threshold is then 



(0 120081) . We adopt a 



aul et alTl201l[). F or JANUS 



( Falcone et al I2OO9I) . The lumi 



47rd?Fii, 



(22) 



Here di is the luminosity distance for the adopted ACDM 
cosmology. EXIST is expected to be 7 — 10 times more 
sensitive than Swift |Grindlavll2010l) . We set the EXIST 
sensitivity threshold to ten times lower than Swift as an 
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approximate estimate. For simplicity, we assume that the 
spectral energy distribution (SED) peaks at X-to-7 ray 
energy (detector bandwidth) as an optimistic case. The 
rate would be less if the SED is very different from Pop 
II/I GRBs. 

2.6. Initial mass function and GRB formation efficiency 

The stellar initial mass function (IMF) is critically impor- 
tant for determining the Pop III GRB rate. We define the 
GRB formation efficiency factor per stellar mass as 



VOBB = /( 



GRB ; ; ' 

Ja/"'' m(j)[m)dm 



(23) 



where (/)(m) is the stellar IMF, and /grb = 0.001 is the 
GRB frac tion, since we expect on e GRB every 1000 su- 
pernovae ( Langer fc NormanI 2006[) . We assume that Pop 
III GRBs have a similar fraction to the standard case. 



Izzard et al.l (|2Q04 ) argue that GRB formation efficiency 
could increase by a factor of 5-7 for low-mctallicity stars 
(^ IO^^Zq). If most of the fi r st sta rs arc rotating rapidly 
as suggested by IStacv et al.l (j201ll) , we can expect that 



a significant fraction of the first stars can produce GRBs. 
Thus given the uncertainty in the parameter for free metal 
stars, we also explore the possibility that /grb is 10— 100 x 
higher. 

We consider the following two forms of IMF. One is a 
power law with the standard Salpeter slope 

<^(m) cx m-2■^^ (24) 

and the other is a Gaussian IM F ( Scannapieco et aDl2003t 
Nakamura fc Umemurall2001 ): 



(m)m ^ dm 



2'K(Jc{M) 



g-(™-M)V2..(A/)^dm. (25) 



For the latter, we assume M = 55OM0 for Pop III.l 
and M = 55M0 for Pop III. 2, with dispersion CTc = 
(M — il/iow)/3. Afiow is the minimum mass for a given stel- 
lar type, IOOMq for Pop III.l, and IOM0 for Pop III.2, 
whereas Mup is the maximum mass for a given stellar 
type, IOOOMq for Pop III.l, and - lOOAf© for Pop III.2. 
-^GRB is the minimum ma ss that is able to trigge r GRBs, 
which we set to be 25M0 (jBromm fc Loebll2006l ). Not all 
Pop III.l stars will leave a black hole behind at their 
deaths. In the narrow mass range of ^ 140 — 26OM0 
Pop III stars are predicted t o undergo a pair-instabil ity su- 
pernova (PISN) explosion (iHeeer fc Wooslevll2002[) . This 
range of mass is excluded from the calculation of Eq. ((23|) . 

The efficiency factor for the power-law (Salpeter) IMF 
is ??grb//grb - 1/926Mq1 and l/STM"! for Pop III.l 
and Pop III. 2, respectively. Using the Gaussian IMF, 
?7grb//grb - 1/538Mq1 and 1/53Mq1 for Pop III.l 
and Pop III. 2. respectively. Thus, the GRB formation ef- 
ficiency for Pop III. 2 can be about an order of magnitude 
higher than Pop III.l because of the lower characteristic 
mass of Pop III. 2 stars. 



3. Redshift distribution of GRBs 

Over a particular time interval, Atobs, in the observer rest 
frame, the number of observed GRBs originating between 
rcdshifts z and z + is 



Q^''gr.b 
dz 



^r,obs 
*GRB 



Atobs dV 



1 



dz ' 



(26) 



where dV/dz is the comoving volume element per unit 
redshift, given by 



dF 
d7 



A-Kcdl 



dt_ 
dl 



(1 + ^) 

Figure [5] shows the intrinsic GRB rate 



diVGRB , .Aiobsdy 

= *grb(^;)- 



dz 



1 + z dz 



(27) 



(28) 



In this plot, we have not considered observational effects 
such as beaming and instrument sensitivity; namely, we 
set riobs = 47r, 77boam = 1, and Lum{z) = in Eq. We 
show the GRB rate for our choice of two different IMFs. 
Interestingly, Fig. [5] shows that the results depend only 
weakly on the choice of IMF. 

Figure [5] shows the most optimistic case, assuming 
a high star formation efficiency /* = 0.1 for Pop III.l, 
/» = 0.01 for Pop III. 2, an inefficient chemical enrich- 
ment, Uwind — 50 km/s, fcRB = 0.1, and a Gaussian IMF 
for both Pop III.l and Pop III. 2 stars. We note that con- 
straints on these quantities will be useful for placing upper 
limits on the GRB observed rate. 

3.1. Radio afterglows 

Follow-up observations of high- redshift GRBs can be done 
by observing their aft e rglows, especially in radio band 
(jloka fc MeszaroEl[200l llnoue et al.ll2007l ). We calculated 
the radio afterglow light curve s for Pop III GRBs follo wing 
the stand a rd pre scription from Sari et al. ( 1998 . 1999() and 
Meszarod (|2006l ). The afterglow light curve at the time td 
is given by the shock radius and the Lorentz factor 7^. 
These two quantities are related by E'iso ~ iTrr^j^nnipC^ 
and Td ^ C7did, where n is the medium density and nip the 
proton mass. The true energy is given by i^truc = ^^£'iso/2, 
where 9 is the half opening angle of the shock. The spec- 
trum consists of power-law segments linked by critical 
break frequencies. These are fa (the self absorption fre- 
quency), fm (the peak of injection frequency), and Vc (the 
cooling frequency), given by 

-coc (l + .)-/^e-/^n-^;^v^t-/^ 



^l,max OC (1 + z)e^J^n^^'^Eisod 



-2 
L ' 



(29) 



where g{p) ~ (p — 2)/(p — 1) is a function of energy spec- 
trum index of electrons (A^(7e)rf7e fx j'^d^je, where 7e 
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is the electron L orentz factor) , and ee and are the ef- 
ficiency factors ( Meszaro^ 2006 ) . There are two types of 
spectra. If < Vc, we call it the slow cooling case. The 
flux at the observer, F^, is given by 

.max 5 ^ ^ ^ai 

{vJvn,)-(P-^y^{vlv,)-P'^F,,^,^, V > V,. 

(30) 

where i^^^max is the observed peak flux at distance cIl from 
the source. 

For t^ni > called the fast cooling case, the spectrum 

is 



F„ 



{vIvcfl^F,^ 
{vIvcY^I-'F,, 



',max; 
'.max? 



max 1 ^ ^ ^m • 



(31) 

As the GRB jet sweeps the interstellar medium, the 
Lorentz factor of the jet is decelerated. When the Lorentz 
factor drops below 6*"^, the jet starts to expand sideways 
and becomes detectable by the off-axis observers. These 
afterglows are not associated with the prompt GRB emis- 
sion. Such orphan afterglows are a natural consequence 
of the existence of GRB's jets. Radio transient sources 
probe the high-energy population of the Universe and 
can provide further constraints on the intrinsic rate of 
GRBs. Even if the prompt emission is highly collimated, 
the Lorentz factor drops 7d < around the time 

1/3 / /) X 8/3 



2.14 



5 X 1054 



0.1 



,-1/3 



(1 + z) days, 



(32) 

and the jet starts to expand sideways. Finally the shock 
velocity becomes nonrelativistic around the time 

^NR- 1.85x10^ (^If^)''' (/j)'''n-/^(l+.)days, 

(33) 

( loka fc Meszaroal2005 ). After time tg, the temporal de- 
pendence of the critical break frequencies should be re- 
placed by Vc oc , Vn\ oc t^^, Us. cx i^i/^, and i^,y,max oc i^i 
(|Sari et al.l[l999l) . We also used the same evolution in the 
nonrelativistic phase for simplicity, which underestimates 
the afterglow flux after ^nr- 

Figure [7] shows the light curves for a typical GRB from 
Pop III. 2 stars assuming an isotropic kinetic energy E{so ^ 
IQ^ipypg ^jjj proportion to the progenitor mass) as a lower 
limit. Pop III.l afterglows are expected to be brighter. 
Consistently with previous works, we conclude that it is 
possible to observe the GRB radio afterglows with ALMA, 
LOFAR, EVLA, and SKA. 



3.2. Upper limits from radio transient survey 

In this section, we derive upper limits on the intrinsic 
GRB rate (including the off-axis GRB) using ~ 1 year 



timescale radio variability surveys. T here are several ra - 
dio transient surveys completed so far. lBower et al.l (|2007| ) 
used 22 years of archival data from VLA to put an upper 
limit of ~ 6 deg~^ for 1-year variability transients above 
90 A tJy, which is equivalen t to < 2.4 x 10^ for the whole 
sky. iGal-Yam et all (|2006l) used FIRStEI and NVS^l ra- 
dio catalogs to place an upper limit of ~ 70 radio orphan 
afterglows above 6 mJy in the 1.4 GHz band over the en- 
tire sky. This suggests less than 3 x 10^ sources above 
0.3 mJy on the sky, because the number of sources is ex- 
pected to be proportional to flux limit f ^J?^^ (assuming 



Eucli dian space and no source evolution) ()Gal-Yam et al.l 
I2OO6I ). From Fig.El a typical GRB's radio afterglow with 
isotropic kinetic energy Eiso ^ lO'^^ ergs stays above 0.3 
mJy over ^ 10^ days. 

By combining the results shown in Figs. [5] and [51 we 
expect ^ 30 — 3 X 10^ sources (10^ — 10^ events per year 
X 10^ days) above ^ 0.3 mJy. (We integrate the event 
rate over redshift.) As a consequence, the most optimistic 
case for Pop III. 2 should already be ruled out marginally 
by the current observations of radio transient sources, if 
their luminosity function follows the one assumed in the 
present paper. Only more conservative models are then vi- 
able. Radio transient surveys are not yet able to set upper 
limits on the Pop III.l GRB rate. The above conclusion 
is model dependent, because the afterglow flux depends 
on the still uncertain quantities, such as the isotropic en- 
ergy £'iso and the ambient density n. If the circumburst 
density is higher than usual, the constraints from the ra- 
dio transient surveys would be even stronger. Also the 
GRB formation efficiency and the beaming factor are not 
known accurately, which can affect both the intrinsic and 
observed rate more than one order of magnitude. 

In Figs. [8]|9l we show the predicted observable GRB 
rate dN^%g/dz in Eq. (gH) for Pop III.l and III.2 de- 
tectable by the Swift, SVOM, JANUS, and EXIST mis- 
sions. The results shown are still within the bounds of 
available upper limits from the radio transient surveys. 
Overall, it is more likely to observe Pop III. 2 GRBs than 
Pop III.l, but the predicted rate strongly depends on the 
IGM metallicity evolution, the star formation efficiency 
and GRB formation efficiency. The dependence on the 
IMF is relatively small. 

Figure [TU] shows the GRB rate expected for EXIST ob- 
servations. Because the power index of the LF is uncertain 
at the bright end, we added two lines to show the resulting 
uncertainty in our prediction. We use the maximum rate, 
which is within the constraints by the current observations 
of radio transients. We expect to observe ^ 20 GRBs 
per year at z > 6 for Pop III. 2 and N ~ 0.08 per year 
for Pop III.l at z > 10 with the future EXIST satehite 
at a maximum. Our optimist case predicts a near-future 
detection of Pop III. 2 GRB by Swift, and the nondetec- 
tion so far could suggest a further upper limit or difference 
between the Pop III and present-day GRB spectrum. 
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Fig. 5. The intrinsic GRB rate dA^^GRB/dz. The number of 
(on-axis + off-axis) GRBs per year on the sky in Eq. ((28|) . 
as a function of rcdshift. We set /, = 0.001, fcBB ~ 0.01 
and Wwind = lOOkm/s for this plot. Salpeter IMF, dashed 
black line, Gaussian IMF, dotted black line, for Pop III. 2; 
and Salpeter IMF, dashed blue line, Gaussian IMF, dotted 
blue line, for Pop III.l. 



105 



Optimistic case 
'\ Pop III.2 



Pop III.l 



Fig. 6. The intrinsic GRB rate dA^GRe/dz. The number 
of (on-axis -I- ofF-axis) GRBs per year on the sky in Eq. 
(P5|) . as a function of rcdshift for our optimistic model. We 
assume a high star formation efficiency; /* — 0.1 for Pop 
III.l; and /* = 0.01 for Pop III. 2; slow chemical enrich- 
ment, Wwind — 50km/s; high GRB formation efficiency, 
fcRB = 0.1; and a Gaussian IMF; for both Pop III. 2, 
dashed black line; and Pop III.l, dotted blue line. 



4. Conclusion and discussion 

There are still no direct observations of Population III 
stars, despite much recent development in theoretical 
studies on the formation of the early generation stars. 
In this paper, we follow a recent suggestion that massive 
Pop III stars could trigger collapsar gamma-ray bursts. 
Observations of such energetic GRBs at very high red- 
shifts will be a unique probe of the high-redshift Universe. 
With a semi-analytical approach we estimated the star 
formation rate for Pop III.l and III. 2 stars including all 
relevant feedback effects: photo-dissociation, reionization, 
and metal enrichment. 

Using radio transient sources we are able to derive con- 
straints on the intrinsic rate of GRBs. We estimated the 
predicted GRB rate for both Pop III.l and Pop III. 2 stars. 



ALMA 
-EVLA 



-ij. 1.4 GHz 

\ \\500 MHz 

' c- LOFAR\ \\ 



\ - - - - 



1 10 



100 
t (days) 



Fig. 7. The theoretical light curve of radio afterglow of 
a typical Pop III. 2 GRB at z ~ 10. We show the evo- 
lution of afterglow flux i^(mJy) as a function of time 
t (days) for typical parameters: isotropic kinetic energy 
Eiso = 10^^ erg, electron spectral index p = 2.5, plasma 
parameters Eo ~ 0.1, ee = 0.01, initial Lorentz factor 
7d = 200, interstellar medium density n = 1 cm""^, for 
the range of frequencies: 500 MHz (dashed brown line), 
1.4 GHz (dashed red line), 10 GHz (dashed black line), in 
comparison with flux sensitivity Fj^™ as a function of inte- 
gration time, tint (days) for SKA (dot-dashed green line), 
EVLA (dot-dashed orange line), LOFAR (dot-dashed blue 
line) and ALMA (dot-dashed purple line). 



Pop III.l 



EXIST 




JANUS 
SVOMf 

0.001 ' 
5xl0-4''/>''Swift 
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Fig. 8. Predicted Pop Ill.lobservcd GRB rate. Those ob- 
served by Swift, dashed red line; SVOM, dot-dashed black 
line; JANUS, dotted blue line; and EXIST, green hue. We 
adopt a GRB rate model that is consistent with the cur- 
rent upper limits from the radio transients; Gaussian IMF, 
t'wind = 50km/s, = 0.1, foRB = 0.1. 



and argued that the latter is more likely to be observed 
with future experiments. We expect to observe maximum 
of A^ < 20 GRBs per year integrated over at z > 6 for Pop 
HI. 2 and A^ < 0.08 per year integrated over at z > 10 for 
Pop III.l with EXIST. 

We also expect a larger number of radio afterglows 
than X-ray prompt emission because the radio afterglow 
is long-lived, for ~ 10^ days above ~ 0.3 mJy from Fig. [71 
Combining with the intrinsic GRB rate and constraints 
from radio transients, we expect roughly ^ 10—10^ 
radio afterglows above > 0.3 mJy already on the sky. 
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Fig. 9. Predicted Pop III.2 observed GRB rate. Those 
observed by Swift, dashed red hne; SVOM, dot-dashed 
black hne; JANUS, dotted blue line; and EXIST, green 
line; for our model with Salpeter IMF, Uwind = lOOkm/s, 
= 0.01, fcBB = 0.01. 




Fig. 10. Predicted maximum GRB rates observed by 
EXIST. We adopted Salpeter IMF, Wwind = lOOkm/s, 
Igrb = 0.01, /* = 0.01; for Pop III.2 GRBs, dashed 
black hne; and Gaussian IMF, Wwind = 50km/s, = 0.1, 
fcRB = 0.1; for Pop III.l GRBs, dotted blue line. Dotted 
red lines represent the same with LP's bright end power 
law index 1.7 and 0.8. 



They are indeed detectable by ALMA, EVLA, LOFAR, 
and SKA, and have even been detected already by ~ 1 
yr-timescale variability surveys. We showed that using a 
semi-analytical approach combined with the current sur- 
veys, such as NVSS and FIRST, we are already able to 
constrain the Pop III. 2 GRB event rate. 

Finally, it is important to note that our knowledge 
of the first stars and GRBs is still limited, and there 
are uncertainties in their properti es, most sign i ficant ly in 
their charac t eristi c mass. Recently, Clark et al.l ( 201lh and 
(j201lh have performed cosmological simula- 



Greif et al. 



tions using a sink particle technique to follow the evo- 
lution of a primordial protostellar accretion disk. They 
find that the disk gravitationally fragments to yield mul- 
tiple protostellar seeds. Although the final mass distri- 
bution of the formed stars is still uncertain, formation 
of multiple systems, especially massive binary Pop III 
stars, would increase the rate of high-redshift GRBs (e.g.. 



iBromm fc L^l200fil: lEV^er fc Hegeijl200,4 and references 
therein). If the GRB fraction per collapse /grb h^ Eq. 
(|23p is much larger than the current one, say /grb ^ 1, 
the Pop III.l GRBs might also become detectable with 
the radio telescopes (~ 300 afterglows above ~ 0.3 mJy 
on the sky) and the X-ray satellites 1 event per year 
for EXIST) in the future. 
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